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ABSTRACT 

The double pulsar system J0737— 3039 is not only a test bed for General Relativity and theories 
of gravity, but also provides a unique laboratory for probing the relativistic winds of neutron stars. 
Recent X-ray observations have revealed a point source at the position of the J0737— 3039 system, 
but have failed to detect pulsations or orbital modulation. Here we report on Chandra X-ray Ob- 
servatory High Resolution Camera observations of the double pulsar. We detect deeply modulated, 
double-peaked X-ray pulses at the period of PSR J0737— 3039A, similar in appearance to the observed 
radio pulses. The pulsed fraction is ~ 70%. Purely non-thermal emission from pulsar A plausibly 
accounts for our observations. However, the X-ray pulse morphology of A, in combination with pre- 
viously reported spectral properties of the X-ray emission, allows the existence of both non-thermal 
magnetospheric emission and a broad sinusoidal thermal emission component from the neutron star 
surface. No pulsations are detected from pulsar B, and there is no evidence for orbital modulation or 
extended nebular structure. The absence of orbital modulation is consistent with theoretical expecta- 
tions of a Poynting-dominated relativistic wind at the termination shock between the magnetosphcre 
of B and the wind from A, and with the small fraction of the energy outflow from A intercepted by 
the termination shock. 

Subject headings: stars: neutron — pulsars: individual (J0737— 3039A, J0737— 3039B) — X-rays: stars 



1. BACKGROUND 

Binary neutron star systems are rare, and even among 
them, the double pulsar system J0737— 3039 is extraordi- 
nary, since both the neutron stars are detected as radio 
pulsars. Th e system consists o f the recycled 22.7 ms 
pulsar "A" (iBureav et all [2 003) and the young 2.8 s 
pulsar "B" (|Lvne et al J 12004 ), in a 2.454 hr eccentric 
(e = 0.09) binary orbit which happens to be nearly edge- 
on to us. As well as being a test bed for General Rela - 
tivity and theories of gravity (e.g. iKramer et alj [2006) . 
the double pulsar is rich in observational phenomena, 
including a short eclipse of A by the magnetosphere of 
B and orbital modu lation of the radi o flux of B due to 
the influence of A (|Lvne et all 12004 ) . The individual 
pulses from B show drifting features due to the impact of 
the low-frequency electromagnetic wave in t he relativis- 
tic wind from A (|McLaughlin et all l2004bl ), while the 
ecl ipse of A is modulated at half the rotational period of 
B (jMcLaughlin et ail [20 04c) . Clearly, the two neutron 
stars have both gravitational and electromagnetic inter- 
actions with each other, and the double pulsar system 
should provide a unique laboratory to investigate the in- 
teractions between the magnetospheres and relativistic 
winds of the two pulsars. 

In this context, the dete ction of X-ray emission 
from the J0737-3039 system (|McLaughlin et alj|2004at 
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iPellizzoni et all 12004 ICampana et all 12004 ) is particu- 
larly exciting. Energetic pulsars generate several forms 
of X-ray emission: quasi-blackbody emission from the 
cooling neutron star surface and/or from heated polar 
caps; pulsed non-thermal emission from the pulsar mag- 
netosphere; and at larger distances from the pulsar, syn- 
chrotron emission from a pulsar wind nebula (PWN) 
powered by the relativistic particle outflow. All of these 
processes may be taking place in the J0737— 3039 sys- 
tem (see, e.g.. iKargaltsev et al.l l2006h . Specifically, the 
X-rays could be pulsed magnetospheric or thermal emis- 
sion from pulsar A (as seen for several other recycled 
pulsars; seelZavlm 2006) , could or iginate in the colliding 
winds of A and B (|Lvutikovll2004D . or could be produced 
by the shock generated when one or both of t he pulsar 
winds interacts with the inters tellar medium (|Lvutikovl 
[200llGranot fc Meszaros1l27)0l . 

The electrodynamics of pulsar winds have been stud- 
ied in considerable detail through the extended PWNc 
typically seen around yo ung and/or high- velocity pulsars 
(jGaensler fc Slanel I2006D . In systems such as the Crab 
Nebula, the PWN is an expanding synchrotron bubble 
centered on the pulsar. Such nebulae act as calorimeters, 
revealing the geometry and energetics of the pressure- 
confined outflow and its termination shock. However, the 
termination shocks seen in such PWNe are typically at 
distances ~ 10 6 — 1Q 9 Rlc from their pulsars (where the 
light cylinder radius of a pulsar rotating at a frequency 
/ is Rlc = c/27r/). In contrast, the two neutron stars in 
the double pulsar system are separated by < 10 3 Rlc,a 
and only 6.6Rlc,b] a termination shock between them 
can thus probe the properties of a pulsar's relativistic 
wind at smaller separations from the central engine than 
ever studied before. Additionally, detection of an or- 
bital phase de pendence in the X-ray emission might be 
expected fe.g.. lArons fc Tavanil [1993T ) . Such variability 
could constrain the geometry of the emission site, thus 
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providing new insights into the wind physics close to the 
pulsar. 

Here we report on Chandra observations of the double 
pulsar which have high enough time resolution to test for 
pulsations from either pulsar and for orbital variability, 
the latter of which might be expected in the bow shock 
or colliding winds interpretation. Forming histograms of 
count rates as a function of phase, we test for modulated 
X-ray emission at the periods of pulsar A and pulsar B, 
as well as for orbital modulation. 

2. OBSERVATIONS AND DATA ANALYSIS 

The J0737— 3039 system was observed with the Chan- 
dra X-ray Observatory using the High Resolution Cam- 
era (HRC-S) in "timing mode" , which provides the high- 
est available time resolution, with events corrected for 
the instrumental wiring error and time-tagged to 16 /is 
accuracy. The observations spanned 10.5 binary or- 
bits but were split into two segments for spacecraft op- 
erational reasons. The first segment of 55 ks began 
on 2006 February 28, while the second segment be- 
gan ~ 67 ks after the end of the first, and spanned 
38 ks. The pulsar system was unambiguously de- 
tected as a point source in both segments, at a position 
07 h 37 m 51 i ;22 -30°39'40"3 (J2000), consistent with po- 
sitions previously determined at X-ray and radio wave- 
lengths ([McLaughlin et alJl2004aHChatteriee et alj|2005t 
iKramer et al.l I2006D at the ~075 pointing accuracy of 
Chandra. 

X-ray photons were extracted from a 1" radius cir- 
cle at the detected position of the J0737— 3039 system, 
and the times of arrival for the photons were corrected 
to the solar system barycenter using the JPL planetary 
ephemeris DE405. Of the 411 photons extracted, we es- 
timate that ^16 counts were contributed by the X-ray 
background. Of course, we cannot identify which of the 
extracted photons came from the background, and nor 
can we assign photons to the individual pu lsars. Instead, 
we use Tempo 7 and timing solutions from IKramer et all 
(2006) to calculate the binary orbital phase and the ro- 
tational phases of both pulsars A and B at which each 
photon was emitted. In keeping with our request that 
the observation be split (if necessary) into integer orbit 
blocks, the orbital phase of the extracted photons ranged 
from 0.840 to 6.854 in the first segment, and from 0.497 
to 4.531 in the second. Thus only ^5% of the orbit was 
sampled 11 times, while the rest was uniformly sampled 
10 times. In the analysis reported below, we have ig- 
nored the minor oversampling, but we have verified that 
discarding 3 detected photons to force uniformity in or- 
bital coverage does not affect our results. 

3. PULSATIONS FROM PSR J0737-3039A 

Forming a histogram of count rate as a function of the 
rotational phase of A, we detect X-ray pulses from pulsar 
A, as illustrated in Figure [TJ The uncertainties on each 
bin (here and elsewhere in this work) ar e la (68%) con - 
fidence intervals estimated according to iGehrelsl (l986). 
The pulsations are double-peaked and deeply modulated, 
with a pulsed fraction / = (Max - Min)/(Max + Min) 
of 0.74^q 1 4. To estimate the significance of the detec- 
tion, we calculate the Pearson x 2 statistic for the pulse 

7 http:/ /www. atnf.csiro.au/research/pulsar/timing/tcmpo 



profile with 16 bins (degrees of freedom i/=15), and find 
X 2 ' jv — 7.05, corresponding to a probability of only 10~ 15 
(~ 8a) that the profile is drawn from a uniform distri- 
bution. 

A visual comparis on of the radio pulse profile of 
PSR J0737-3039A ([Manchester et al.l 120051 ) with the 
X-ray profile shows a distinct resemblance (Figure [I}. 
iDemorest et al.l (|2004f ) model the radio pulse as two cuts 
through a wide cone of emission centered on a single 
magnetic pole of A, which has its spin and magnetic 
axes nearly aligned (4° ±3°). Although a wide range 
of misal ignment is currently pe rmitted by radio obser- 
vations (|Manchester et al. 2005), both peaks of the ob- 
served pulse appear to come from one magnetic pole, 
implying a very wide fan beam in some geometries. The 
X-ray pulse profile also shows two peaks, whose locations 
fall within the range in pulse phase delimited by the radio 
peaks when the pulses are phase- aligned. 8 This suggests 
that the X-ray emission is from a narrower cone than 
the radio beam. Specifically, the peaks in the X-ray pro- 
file (Figure [1]) are located at pulse phases <j) A ~ 0.27 
and ~ 0.77, as estimated by binning the observed X- 
ray photons at various resolutions, and the peak-to-peak 
separation is - 0.50 ± 0.01 (182° ± 3°), while the p eaks 
in the 1.4 GHz radio profile (Manche ster et aill2005f ) are 
at <j) A = 0.234 and <f> A = 0.789, separated by ~ 200°, 
and the outer rims of the radio emission profile are at 
4>A = 0.164, 0.836. The X-ray emission also shows a sig- 
nificant "bridge" between the two peaks, implying that 
the cone of X-ray emission is (partially) center-filled in 
this model, unlike the broader, hollow radio emission 
cone. 

The detected pulses are quite unlike the typical X- 
ray emission observed from other recycled pulsars with 
comparable spin parameters (e.g. , PSR J0437— 4715 , 
iZavlin et all l2002t iBogdanov etatl 120061 : iZavlinl I2006T) . 
which show broad, roughly sinusoidal pulsations with 
low pulsed fractions and thermal spectra. Instead, the 
pulsations from pulsar A resemble non-thermal pulses 
seen only from t he most energetic recy cled pulsars (e.g., 
PSR B1821-24. iRutledge et all [200l . even though A 
is slower rotating and has a lower spindown energy loss 
rate (E A — 5.9 x 10 33 erg s" 1 ). In this context, we 
note that the fo rmation scen ario for double neutron star 
binary systems ()Stairsl[2004h can result in a shorter ac- 
cretion episode and thus a higher surface magnetic field 
strength compared to other recycled pulsars. Pulsar A 
has an inferred dipole magnetic field B A = 6.4 x 10 9 G, 
comparable to PSR B1821— 24 and the binary pulsar 
B1534+12, but significantly higher than other recycled 
pulsars that show predominantly thermal X-ray emis- 
sion. Both pulsars B1821-24 and J0737-3039A also 
lie above the death line for curvature radiation esti- 
mated bv lHarding et al.1 (|2005l ). suggesting that the pro- 
cesses that power non-thermal magnetospheric emission 
in PSR B1821— 24 may also operate for pulsar A, al- 
though the two differ substantially in period and E. 

The absence of any useful energy resolution in Chandra 

8 IRut lcdgc et al. ( 20Q4|) show that absolute phase alignment is 
possible at the 60 fis level between HRC-S and radio observations 
of the recycled pulsar PSR B1821— 24. Since we have to predict 
and account for both orbital and rotational phase, our timing errors 
are somewhat larger, but insignificant compared to the bin width 
of ~1.4 ms. 
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Fig. 1. — X-ray pulse profile of pulsar A, obtained by folding 89 ks of Chandra HRC-S data with Tempo and a DE405 timing s o lution . 
The uncertainties on each bin (here and elsewhere in this work) are la (68%) confidence int ervals estimated accord ing to Gchrcls (1986). 
The pulse profile is shown twice for clarity, and a radio pulse profile obtained at 1.4 GHz (Manchester et al. 2005) is plotted below (in 
arbitrary units) for comparison. Both profiles are folded using the radio timing solution, and are therefore aligned in phase. 



HRC data precludes spectral fits to the data, but pre- 
vious Chandra ACIS observations can be well-modeled 
by a power law with a photon index T ~ 2.9 ± 0.4 
(|McLaughlin et alj l2004al) . and XMM data is well-fit 
by a power law with a photon index T ~ 3.5^Q 3 
(|Pellizzoni et alj I2004D. Joint fits to the Chandra and 
XMM data (|Campana et all 120041 ) allow for both power 
law (r = 4.2±2-i) and thermal black body (kT bb = 0.20 ± 
0.02 keV) interpretations. Similar fit p arameters (T ~ 3 
or kT bb ~ 0.2 keV) were fou nd by iKargaltsev et all 
:2(l(l(); as well. Additionally, ICampana et all (|2004D 
show that a two-component fit with a fixed power law 
index T — 2 and a black body component (kTbb — 
0.16 ± 0.04 keV) is consistent with the Chandra ACIS 
and XMM data, although two components are not sta- 
tistically required. 

The X-ray spectrum, in combination with our detec- 
tion of sharp, double-peaked X-ray pulses, is thus con- 
sistent with a purely magnetospheric origin for the X- 
ray emission, but it is also possible that the observed 
X-ray pulsations consist of both non-thermal magne- 



tospheric emission and broad sinusoidal thermal pul- 
sations from the hot polar cap. The pulse profile of 
A shows a floor of X-ray emission (Figure [T|), corre- 
sponding to a count rate of w 1.5 ± 0.6 cts ks -1 at 
every phase. An image of the off-pulse counts reveals 
no extended nebular structure, and their distribution 
is consistent with the on-pulse photons. Other recy- 
cled pulsars also show emission at all pulse phases, 
whether their puls ations are broad a n d thermal (e.g., 
PSR J0437-4715 JZavlin et all 120011 iBogdanov etall 
[200l iZavlinl l2006h or narrower and primarily non- 
thermal (e.g., PSR J0218+4232. iKuiper et all 120021 : 
PSR B1821-24. fRutledge et all l2004h. Such unpulsed 
emission is usually ascribed to thermal X-rays emitted 
from the neutron star surface. Assuming that the entire 
X-ray flux of the double pulsar system arises only from 
the combined thermal and non-thermal emission from 
PSR J0737— 3039A, we find that the maximum ampli- 
tude sinusoid A(l + sin27r(0 — 4>o)) that is consistent 
with the observed profile at la could account for as 
much as ~ 60% of the observed X-ray counts, although 
the actual fraction is likely to be far lower. Rotational 
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phase-resolved spectroscopy with substantially more X- 
ray counts will be required to verify or rule out such a 
two-component model. 

In order to investigate possible orbital variations in 
the X-ray pulse profile of A, 9-bin pulse profiles were 
constructed for each quadrant of the orbit. Each of the 
four profiles was then compared to the pulse profile con- 
structed by averaging the other three quadrants. The 
resulting x 2 jv values range between 0.8 and 1.3 (with 
v — 9 degrees of freedom), consistent with no variations. 
While we lack the S/N to definitively rule out any differ- 
ences between the X-ray pulse profiles, no orbital varia- 
tions axe_detected injhe radio pulse profiles of A either 
(e.g. lKramer et al.ll2006l ). 

We note in passing that our estimate of the pulsed frac- 
tion / = 0.74^q'^ is only marginally consistent with the 
upper limit of 60% on the pulsed fraction (assum ing si- 
nusoidal pulses) inferred by iPellizzoni et all (|2004l ) from 
XMM-Newton observations. Since the detected pulse is 
non-sinusoidal, a direct comparison is not possible, but 
~ 60% of our detected photons are above the estimated 
minimum count rate baseline, and > 51% are > 1 a 
above the basel i ne lev el. The XMM pn observations of 
IPellizzoni et al.l (|2004l ). which were in continuous clock- 
ing mode, were totally dominated by the background due 
to the one-dimensional readout, while the XMM MOS 
chips lack the time resolution to detect pulses from A, 
leading to a limit which is less robust compared to the 
Chandra HRC detection presented here. 

4. NON-DETECTION OF PSR J0737-3039B 

We repeated the analysis described in § 3 for PSR 
J0737-3039B. The results are shown in Figure d No 
X-ray pulsations are detected, either by folding the full 
span of data, or by selecting counts which are in the 
lowest emission bins of the pulse profile of A, 0.46875 < 
4>A < 0.65625. For the folded profile from the entire 
data span, we calculate x 2 l v = 1.34, corresponding to 
a 17% probability that the data ar e drawn from a un i- 
form distribution. As described by iLeahv et al.l (|1983f) . 
epoch folding is not as sensitive to broad, smooth pulses 
as the family of Rayleigh statistics Z^, which also avoid 
the need to bin data. Therefore, we also calculated the 
H statistic (H = M ax(Z 2 , - 4m + 4), for 1 < m < 20; 
Ide Jager et al.l fl989h . which is well suited to searching 
for an unknown modulation shape. We find H = 0.035, 
at m = 1, corresponding to a null hypothesis probabil- 
ity (i.e., the probability that we are sampling a uniform 
distribution) close to unity. 

Pulsar B shows significant enha ncements in radio emis- 
sion at some parts of its orbit (|Lvne et al.ll2004T) . but 
folding X-ray photons selected from those orbital phase 
ranges does not show any evidence for pulsations cither. 
The non-detection is unsurprising, since pulsar B has 
rotational parameters and a spindown energy loss rate 
(Eb = 1-7 x 10 30 erg s^ 1 ) similar to other "ordinary" 
pulsars (ages ~ 10 6 — 10 8 yr), which are not known for 
their X-ray emission. The spindown luminosity of pulsar 
B is only ~3x 10~ 4 Ea, and so pulsar A is expected to 
dominate any X-ray emission from the system. 

5. LIMITS ON ORBITAL MODULATION 

Sys t ems such as the Crab n ebula (Ke nnel fc Coronitl 
119841: iGallant fe Aronsl I1994D and the nebula around 
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Fig. 2.— Non-detection of pulsations from PSR J0737-3039B. 
89 ks of Chandra HRC-S data were folded at the rotational phase 
of B, as predicted by Tempo, but no pulsations were detected. 
Two pulse periods are shown for clarity. Further, we extracted 26 
photons detected in the off-pulse phase range of A, 0.46875 < 4>a < 
0.65625, corresponding to the three bins with lowest photon counts 
in Figure [T] The folded counts are shown on the same scale (the 
lower pulse profile in the figure). Again, no significant pulsations 
were detected. 



PSR B1509-58 (jGaensler et all [20021 ) provide the best 
current constraints on the behavior of pulsar winds at 
large distances (~ 10 6 — 10 9 Rlc) from the neutron 
star, but few constraints exist for the close-in behav- 
ior. The magnetization parameter cr, the ratio of Poynt- 
ing flux to the kinetic energy flux in the wind, is a 
key descriptor of such systems. Optical, near-infrared, 
and X-ray images at sub-arcsecond resolution reveal 
that the shock has an axisymmetric structure of equa- 
torial arcs (wisps ) and polar jets (knots) that vary on 
short time-scales (iHester et aL|[200a iPavlov et al.l 120031: 
iMelatos et al.ll2005[ L and that the wind transforms from 
a Poynting-dominated outflow (cr 3> 1) near the pul- 
sar to a kinetic-energy- dominated outflow (a < 1) at 
the termination shock ([Kennel fe Coronitll Il984) . Re- 
cent work has begun to elucidate the collimation mech- 
anism that produces the axisy mmetric structure (e.g., 
iKomissarov fe Lyub arskv 2004) . while the conversion of 
Poynting flux to mechanical energy remains poorly un- 
derstood. 

The wind interaction of a neutron star with a stellar bi- 
nary companion allows constraints on the wind behavior 
at ~ 10 4 i?LCi an d such interaction can produce radio and 
high energy emission signatures. For example, the Be 
star — pul sar binary B1259 — 63 produces unpulsed radio 
emission (Bal l et aLlll999l) as well as unpulsed high en- 
ergy emission (e.g.. ICrrove et al.lll995f) . which arise from 
the shock formed between the ste llar outflow and the pul- 
sar wind ([Tavani fc Aronsl [T997). The interaction of the 
pulsar B1957+20 with its white dwarf binary companion 
is expected to pro duce orbital modulation in the X-ra y 
emission (see, e.g.. lArons fc Tavanil[l993l : lMiche]lll99l . 
although the observational evidence for such modulation 
(jStappers et alj 120031: iHuang fc Becker! I2007D is not sig- 
nificant. 
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As opposed to the interaction between a neutron star 
relativistic wind and the particle wind of a stellar com- 
panion, the double pulsar presents a situation where the 
relativistic wind interacts with the magnetosphere of an- 
other neutron star. Additionally, the system separation 
is < 10 3 Rlc,a and only 6.6Rlc,b- The detection of 
orbital modulation in the system would thus be of par- 
ticular interest, since it probes the beh avior of the pulsar 
wind in a high-cr regime. We note that iKargaltsev et al.l 
(2006) find possible orbital phase dependence for the 
double neutron star binary B1534+12, but not for the 
J0737— 3039 system. Their result is based on earlier 
XMM and Chandra ACIS data that lacked the time res- 
olution to detect pulsations from pulsar A. 

Given our detection of deeply modulated pulsed emis- 
sion from PSR J0737— 3039A, we attempted to detect 
orbital modulation in the X-ray emission by folding X- 
ray photons from the off-pulse phase of A, 0.46875 < 
4>A < 0.65625, corresponding to the three bins with low- 
est photon counts in Figure Q] The 25 counts thus se- 
lected from the observation (corresponding to a reduced 
effective exposure of 16.7 ks) were folded as a function of 
orbital phase, and the results are shown in the top panel 
of Figure [31 We find an apparent enhancement at a phase 
~0.69. At the epoch of observation, that phase bin in- 
cludes an orbital longitude u> = 0°, corresponding to A's 
crossing of the ascending node of the orbit. However, 
there is no obvious physical mechanism that could pro- 
duce such an enhancement, and the binned distribution 
has x 2 /v = 1.45, corresponding to a chance probability 
of 11.5%. As in §4, we calculate the H statistic with the 
un-binned orbital phase values. We find H = 7.35 at the 
fifth harmonic Zf , which allows the null hypothesis that 
we have sampled uniformly distributed data at ~ 5%, a 
probability that is small but not insignificant. 

We also check for the enhancement by folding all the 
X-ray photons, and by selecting and folding photons from 
the bridge of emission between the two peaks of A's pulse 
profile (-0.15625 < (f> A < 0.03125), where the contribu- 
tion of the pulsar itself is reduced. The results are shown 
in the middle and bottom panels of Figure [3) and in each 
case, we again calculate the H statistic. H = 1.66 at 
m = 1 when including all the extracted photons, corre- 
sponding to a null hypothesis probability of 52%. For the 
photons chosen from between the two peaks of A's pulse, 
H = 0.023 at m = 1, which allows the null hypothesis at 
a probability close to unity. Together, these results lend 
weight to the conclusion that the apparent orbital mod- 
ulation seen above (with a chance probability of 5%) is, 
in fact, not real. We have also confirmed the absence of 
significant modulation by binning as a function of orbital 
longitude rather than phase, with very similar results. As 
outlined in § 3, it is more likely that the impulsed X-rays 
have their origin in thermal emission from the surface of 
pulsar A. We thus confi rm the negative result found by 
IKargaltsev et all (|2006h . 

Fro m the drifting sub-pulses d etected in B's radio emis- 
sion ([McLaug hlin et aT1l2004bl ). it is apparent that the 
low-frequency electromagnetic wave in the relativistic 
wind from A influences the emission of pulsar B, and 
several models have been proposed where the formation 
of a shock between the two pulsars sho uld produce or- 
bital modulation in their emission (e.g. lLvutikov| [200l 
iGranot fc Meszaros!l2004l : iTurolla fc Trevesll2004f 
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Fig. 3. — Searching for orbital modulation in X-ray emission 
from the J0737— 3039 system. In all cases, the orbit is plotted 
twice for clarity. Top: We extract 25 photons detected in the off- 
pulse phase range of A, 0.46875 < ij>A < 0.65625, corresponding 
to the three bins with lowest photon counts in Figure [T] Folding 
at the binary phase shows an enhancement at an orbital phase 
~0.69. At the epoch of observation, that phase bin encompasses 
an orbital longitude u> = 0°, when A crosses the ascending node 
of the orbit. Middle: Folding all detected photons at the binary 
phase does not show such an enhancement. Note that we lack 
enough counts to detect or constrain the eclipse of A at an orbital 
longitude uj = 90°. Bottom: We extract 65 photons in the mid- 
pulse of A, —0.15625 < <f>A < 0.03125, corresponding to the three 
bins between the peaks of the profile in Figure [T] The absence of 
any significant enhancement confirms that the apparent signal in 
the top panel is spurious. 



However, only a small fraction of the wind power emit- 
ted by A (and half of the power emitted by B) is in- 
tercepted by the shock between the two pulsars, re- 
ducing proportionately the maximum X-ray flux that 
the shock emits. For example, if we assume that the 
wind energy is radiated isotropically from A, and that 
it is intercepted by a sphere centered on B with ra- 
dius Rlc,b, then the power intercepted by the shock, 
E s = 0.006£U + 0-5E B ~ 0.006.E A . If, instead, A's wind 
is intercepted at the surface where pressure balance is 
achieved between the wind from A and the magneto - 
sphere of B, at ~ 0.20 lt-s from B (jLvne et al.ll2004h . 
then we have E s = O.OOI^a + 0.5^ B « O.OOli^- Fi- 
nally, if the shock roughly coincides with the region cen- 
tered on B that eclipses the radio pulses from A, we have 
E s = lQ- b E A + 0.5E B « 1.5 x lQ- i E A (although the 
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processes that contribute to radio eclipses are likely to 
be quite different from those that cause X-ray emission) . 

Of course, the wind radiated from A is unlikely to 
be isotropic, especially if the magnetic and rotational 
axes are nearly-aligned (jDemorest et al.l l2004f ) . and the 
shock geometry is not described simply by intersecting 
spheres centered on A and B. Nevertheless, the conser- 
vative geometric estimates above demonstrate that the 
X-ray power output from the shock is E s < 0.006Ea, 
possibly modulated at the orbital period. Interestingly, 
spectral fits to the Chandra and XMM data imply an X- 
ray efficiency L X /E A < 2 x 1(T 4 (jCampana et al.ll2004T ). 
where L x is the X-ray luminosity in the 0.5 — 10 keV 
range. Thus, if the entire E s were converted to X-ray 
emission, at least two of our proposed scenarios above 
would have resulted in a higher X-ray efficiency for the 
J0737— 3039 system than actually observed. Since we de- 
tect X-ray pulses from A which account for a significant 
proportion (and arguably ~100%) of the observed X-ray 
emission, all of E s evidently does not appear as X-ray 
emission. (We note that for Ea = 5.9 x 1 33 erg s _1 , the 
relations derived bv lPossenti et a l. (2002) for X-ray lumi- 
nosity in the 2 — 10 keV range predict a maximum X-ray 
efficiency L x /Ea < 0.005, consistent with observations.) 

The wind interaction in the double pulsar system 
is fundamentally different energetically from wind con- 
finement in a Crab-like pulsar wind nebula, since the 
termination shock of the wind is much closer to pul- 
sar A (< 10 3 Rlc,a) than in Crab-like nebulae (~ 
10 8 Rlc)- All modern wind models, whether for a 
steady-state, force-free, magne tohydrodynamic outflow 
(jConto poulos & Kazanas 2 0021) or a wave-like, striped 
outflo w (jMelatos fc Melrosd Il996t iLvubarskv fc Kirkj 
|2001[ ). predict values of the magnetization parameter 
a 1 (probably > 100) at these distances, unlike ter- 
mination shocks in pulsar wind nebulae , whe re a -C 1. 
For a high-er shock. iKennel &: Coron iti (1984) estimate 
an upper limit on the power fed into the accelerated elec- 
trons (and hence on the X-ray luminosity of the shock) of 
E s /(8y/a). In summary, as a result of the high expected 
value of a and the small solid angle over which the wind 
from A is intercepted by B, the shock produced at the 
interaction region is unlikely to show significant X-ray 
emission. Similar arguments apply to t he absence of un- 
pulse d radio emission from the system (|Chatteriee et al.l 
l2005h as well. 

6. CONCLUSIONS 



With 89 ks of Chandra HRC observations, we 
have detected deeply modulated emission from 
PSR J0737-3039A. The off-pulse emission reveals 
no extended structure. No pulsations were detected 
from PSR J0737— 3039B, and no orbital modulation was 
detected either. Although we cannot absolutely rule out 
orbital modulation or emission from other mechanisms 
such as bow shocks, we have shown that the entire X-ray 
emission from the J0737— 3039 system can be explained 
as arising from pulsar A alone, either as non-thermal 
magnetospheric emission, or as a combination of mag- 
netospheric and thermal emission. Pulse phase-resolved 
spectroscopy will allow discrimination between these 
two scenarios. 

Like the dog that did not bark in the night, the ab- 
sence of significant and detectable orbital modulation 
in the X-ray emission from the J0737— 3039 system is 
noteworthy. The wind from pulsar A impinges on and 
compresses the magnetosphere of B, leading to deep or- 
bital modulation in the detected radio pulsations from B 
(|Lvne et al.ll2004t ), and the impact of the low- frequency 
electromagnetic wave in the relativistic wind from A 
is also seen in the drifting sub-pulses of emission ob- 
served from B ()McLaughlin et all l2004bf ) . Given the 
strong influence of A on the radio emission from B, it 
may seem natural to ascribe the X-ray emission from the 
PSR J0737— 3039 system to a particle shock formed at 
the wind-magnetosphere interaction site. As we show 
here, such an interpretation is neither favored by theory, 
nor required by the X-ray observations. Our observations 
reveal no significant orbitally modulated shock emission, 
consistent with models for relativistic winds that require 
a Poynting-dominated wind close to the pulsar, and with 
only a small fraction of the energy outflow from A inter- 
acting with the termination shock. 
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